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We compute the time evolution of correlation functions after quantum quenches in the sine- 
Gordon model within the semiclassical approximation which is expected to yield accurate results 
for small quenches. We demonstrate this by reproducing results of a recent form factor calculation 
of the relaxation of expectation values. Extending these results, we find that the expectation values 
of most vertex operators do not decay to zero. We give analytic expressions for the relaxation of 
dynamic correlation functions, and we show that they have diffusive behavior for large timelike 
separation. 


I. INTRODUCTION 


Questions of the relaxation and thermalization of iso¬ 
lated quantum systems hav e attracted a lot of atten¬ 
tion over the last decade 11 ^. Under which conditions a 
given system relaxes or thermalizes? If the asymptotic 
stationary state is not thermal, can it be described within 
the framework of statistical physics? How quickly is the 
asymptotic state reached and what are the characteris¬ 
tics of the time evolution? The increased interest in these 
fundamental questions is to a major part due to the spec¬ 
tacular advances in cold atom experiments which are able 
to study the coherent evolution of isolated quantum sys¬ 
tems, moreover, many properties and parameters of the 
systems are tunable^®. 

With the scope of understanding thermalization or the 
lack thereof, many studies focused on the asymptotic 
steady state after a quantum quenctf-. Still, there are 
many unsettled questions. An example is the General¬ 
ized Gibbs Ensembl^ that was proposed to describe the 
stationary state of integrable systems which possess an 
excessive number of conserved quantities. Lately, its ap¬ 
plicability for continuum systems was questioned®} and 
its natural implementation using local conserved quanti¬ 
ties was even shown to fail to capture the steady state of 
the XXZ spin chaiiPSE] 

Even less is known about the details of the relaxation 
process and the time scales of the relaxation. Numeri¬ 
cal approaches are usually constrained either by the sys¬ 
tem size or the times until they are able to follow re¬ 
liably the evolution of the system, and the long time 
behavior in the thermodynamic limit is very difficult to 
study. This is especially true for continuum many-body 
systems. A promising direction is based on the quench 
action methocP^HHI 

(see below). Progress in the analytic 
description has been made in conformal field theoriesP 
and in s ystem s that can be mapped to free bosons or free 
fermion,} l:j l 15 l EH. A notable exception is Ref. [521 . 

The relaxation of quantum many body systems can 
happen in a number of steps. For example, a weakly non- 
integrable integrable system can reach a prethermaliza- 


tion platearP- 31 i22H26] close to the steady state of the in¬ 
tegrable model and real thermalization takes place only 
at much larger timescales. 

Obtaining analytical results for genuinely interacting 
systems is notoriously hard even in integrable models, 
where the spectrum and the matrix elements of local op¬ 
erators, the so-called form factors, are usually known. 
Based on these ingredients, a linked cluster expansion 
can be constructed where the small expansion param¬ 
eter is esse ntially the density of excitations after the 
quendJ2ZI2SI Summing up the series is a daunting task, 
nevertheless, some progress can be made at least regard¬ 
ing the short time behavior after a quantum quench. 
Another possible approach is based on the so-called 
quench action technique for Bethe Ansatz integrable 
system^ 2 which works directly in the thermodynamic 
limit. Both methods have been applied to study the re¬ 
laxation dynamics in the sine-Gordon model in Ref. [35J 
and predicted an exponential decay of the vertex opera¬ 
tor (e^^’ 4 )/ 2 ) ~ e -4 / T . 

Here we shall follow a complementary and more in¬ 
tuitive approach, and study the evolution of dynami¬ 
cal correlations after a quantum quen ch bv extending 
the semiclassical approach of Refs ] 18 l 22 llHE§]E3 to study 
small quenches in the gapped phase of the sine-Gordon 
model. The sine-Gordon model is a paradigmatic model 
providing the low energy effective description of a wide 
range of one-dimensional systems including spin chains, 
spin ladders, and cold atomic gaseiP^^I. It is defined by 
the action 


4 ( a t $) 2 -(^$) 2 + Acos (/ 3 $) , 
c z 

(1) 

with <l>(x,t) a bosonic field. In this work, for the sake 
of simplicity, we shall focus on the gapped repulsive 
phase (1/ -s/2 < j3 < 1 ), where the model contains mas¬ 
sive topological excitations (kinks), the so-called solitons 
and antisolitons with charge m = ±1 interpolating be¬ 
tween neighboring minima of the cos(/3$) potential, but 
no bound states (breathers) exist 

The time evolution of the vertex operator was 
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studied after quenches in the attractive regime in Ref. [52, 
Analytic results have been derived for the correlations 
of the same operator for quenches between the exactly 
solvable points A = 0 and /3 = 1/ y/2 in Ref. [213 In Ref. 
1541 the time evolution of an inhomogeneous initial state 
was studied. 

The semiclassical method is based on the observation 
that at small temperatures or after a small quench the 
density of quasiparticles as well as their velocity is small. 
Therefore quasiparticles can be treated as classical enti¬ 
ties apart from collisions, where their de Broglie wave¬ 
length becomes comparable to their separation. These 
quantum effects are taken into account by using the low 
energy limit of the two-particle scattering matrix. 

The semiclassical approach turned out to yield a 
remarkably accurate description of many gapped sys¬ 
tems. It has been successfully applied to compute fi¬ 
nite temperature correlation functiond^SHIl 

as well as 

the time evolution of correlations after global and local 
quenche d — I I M 125III1 

As we demonstrate here, this approach is able to cap¬ 
ture the leading behavior of the decay processes and re¬ 
produces the tediously obtained results of Ref. [3B3 with 
ease. 

However, we can go significantly beyond these results. 
While extending the form factor based calculations of 
Ref. [28] to other observables (e.g. two-point functions) 
seems to be a very demanding task, in the semiclassical 
method it poses only slight, surmountable complications. 
Assuming that only soliton-antisoliton pairs are present 
in the initial state (as is the case for fixed d> initial con¬ 
ditions), we obtain new results in two directions. On the 
one hand, we calculate the relaxation of general vertex 
operators 

G a (t) = = (iJ; 0 \e iHt e ia ®( x, °^ e~ iHt \il>o) (2) 

with the somewhat surprising result that — in the uni¬ 
versal limit - they do not decay to zero but approach 
finite asymptotic values. 55 On the other hand, we com¬ 
pute the time evolution of dynamical two-point functions 
of general vertex operators 

C a (x' - X-,t,t') = . ( 3 ) 

We show that the two-time correlations show diffusive 
behavior for generic values of a. This is expected to some 
extent given that diffusive behavior was also observed 
in the semiclassical treatment of correlations in thermal 
equilibriuirP®H. 

The paper is organized as follows. In Section II we de¬ 
scribe the semiclassical approach in detail. Expectation 
values of vertex operators are computed in Section III. 
We calculate the time evolution of general dynamic cor¬ 
relation functions in Section IV, and we analyze the equal 
time correlations, the local correlation functions and the 
correlations in the asymptotic state separately. We give 
our conclusions in Section V. 


m‘ 2 m'-L 



FIG. 1: A two-body collision of kinks. 


II. THE SEMICLASSICAL METHOD 


As discussed in the Introduction, for small quenches, 
the energy density injected in the system is small, and 
quasiparticles are generated with a low density and with 
energies close to the energy gap. In this low density limit, 
the motion of quasiparicles is ‘slow’ and can be treated 
semiclassically. The quantum expectation values ([2]) and 
([3]) are calculated as averages over the kink configura¬ 
tions, that is over the initial positions, velocities and 
charges of the kinkJ- 59 ! Importantly, due to momentum 
and energy conservation, the trajectories of quasiparti¬ 
cles remain straight lines in ID and follow “rays”, (see 
Figs. 02 ). 

In the particular case of the repulsive sine-Gordon 
model, in the semiclassical limit, most of the time the 
field <f> remains close to minima of the cosine term in 
Eq. 0 


<f> = n 


2tt 

J ’ 


n £ Z. 


(4) 


Quasiparticles are just kinks (domain walls) separating 
domains of constant $, such that n increases (decreases) 
by one when crossing in the positive spatial direction a 
soliton (antisoliton) trajectory. 

In the small density limit only two-particle collisions 
are relevant. Given the small velocity of quasiparticles, 
the scattering matrix of solitons and antisolitons can be 
approximated by its low momentum limit as 


= (5) 

i.e. kinks scatter as “hard balls” (c.f. Fig. 0. This struc¬ 
ture is crucial for the rest of this work: it implies that the 
spatial sequence of the topological charges of the kinks 
(solitons and antisolitons) remain the same for all times 
in this asymptotic limit. Alternatively, in terms of do¬ 
mains, Eq. 0 implies that the ‘color’ sequence of do¬ 
mains remains invariant under time evolution (see Fig.0. 

As discussed above, any quantum quench of the Hamil¬ 
tonian generates a gas of quasiparticles. However, very 
importantly, in many cases the post-quench state seems 
to possess a particular, correlated structure in terms of 
quasiparticle pairs. Following Ref. 1381 we take our initial 
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state to be 


{ f°° d 9 j 

J Q ^ K mm'(O)Zl(-0)ZU0)j |0>, (6) 

where Z^O) creates a kink of type m = ± with relativis¬ 
tic rapidity 9. These states are coherent superpositions of 
kink pairs and are called integrable initial states due to 
their resemblance to integrable boundary stated^. The 
exponential form featuring a single creation amplitude 
(for each pair type) generates independent pairs of par¬ 
ticles with opposite velocities. For small quenches, the 
pair creation amplitudes K m m'{9) and the correspond¬ 
ing densities are small, 

c\0 

Pmm' ~ J — Mccosh(e)\K mm '(e)\ 2 < 1 , ( 7 ) 

where M is the kink masd^ . 

For small quenches, the form ([g] of the initial state 
can be justified on very general grounds: by momentum 
conservation, a local perturbation, to the lowest order, 
gives rise to pairs of kinks flying away from each other 
with the same velocity. Moreover, if the field is originally 
constant (Dirichlet boundary condition) then, since the 
field must remain unaltered away from the perturbation, 
the total topological charge of each pair must be zero, i.e. 
pairs must form soliton-antisoliton pairs. Since a homo¬ 
geneous global quench by the integral of a local operator 
is a sum of such local quenches, thus the post-quench 
state will be populated by independent pairs. Indeed, 
this pattern of the quench creating pairs of quasiparti¬ 
cles has been obs erved for several integrable systems and 
initial stateiP^EZHSI], even for large quenches. The variety 
of quenches featuring the pair structure suggests that 
this may be a general phenomenon. We note that the 
sine-Godon model with the initial state © appears in 
the description of interference patterns between split ID 
condensat es an d the Ramsey sequence of two-component 
ID boson^Ml 

We thus assume that the post-quench semiclassical 
configuration consists of a collection of uniformly and 
independently distributed pairs of straight lines, placed 
along the t = 0 axis, and a random sequence of soliton- 
antisoliton charges (see Fig. [2]). Kink pairs have a velocity 
distribution f m m'{v ) (v > 0), i.e. each pair of kinks with 
charges m and m' traveling with velocities —v and v is 
created with a probability density f rnm > (v). For the spe¬ 
cific initial state (16]) this distribution is related to the 
amplitudes K mrn ’(9) as 

fmm'(v) ~ ^-\Kmm'{v/c)\ 2 , (8) 

27 Tp 

where p = m' Pmm' is the total density of pairs 
and we used that the velocities are nonrelativistic and 
v/c = tanh(9 « 9. With this definition, the distribu¬ 
tion functions are normalized as J n dvf mm '(v) = 

mm' Pmm 1 = 1- 



FIG. 2: A kink configuration for initial states with 
soliton-antisoliton pairs only. In this example 
n L = 1 ,n R = l,n 7 = 1 ,n A = 1. 

In our case, only soliton-antisoliton and antisoliton- 
soliton pairs are created with equal probability, thus 
K ++ = K _= 0 and K + _ = K _= K. We can there¬ 

fore characterize the velocity distribution of the kinks by 
a single function, 

f(v)=-\K(v/c)\ 2 , / dvf(v) = 1. (9) 

np Jo 

Our goal is to determine the time evolution of the ex¬ 
pectation value ([2]) and the correlation function ^ . Cor¬ 
relation functions depend on the difference <t>(£i, t\) — 
$(£ 2 , £ 2 ) which, apart from a prefactor 2tt/(3 is equal to 
the sum of charges of the kinks that cross the straight 
line connecting the two points where the operators are 
inserted. As noticed earlier, due to the special structure 
of the universal S-matrix ]5]) , the spatial sequence of do¬ 
mains of constant $ remains unchanged under the time 
evolution (see Fig. [ 2 ]). Domains can thus be labeled by 
an integer, and &{x,t) = in case the point (x,t) lies 
in the ( th domain. Consequently, <F(£i,ti) — $(£ 2 ,^ 2 ) = 
if (£ 1 , ti) and (£ 2 ,^ 2 ) lie in the Z/ h and l^ 1 do¬ 
main, respectively, and 

OjT j 2 

= -j mi ' ( 10 ) 

p i=h +1 

where {m,} are the charges of the s = |Z 2 — (i| domain 
walls (kinks) encountered while going from domain // to 
( 2 - The precise value of s is determined by the number and 
directions of the rays that intersect the segment between 
the two operator insertion points, and an averaging over 
the initial positions {£ 7 } and velocities {///} of the pairs 
as well as over the charges {mj needs to be performed. 
The correlation function ([ 3 ]) is thus expressed as 

C7 Q (A£; t ,0=4° ) (A£,t'-t)^^ 

m, , 

with Ca' 1 denoting the vacuum correlation function. 
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III. RELAXATION OF EXPECTATION VALUES 

Let us start by calculating the time evolution of the 
expectation value . For this we need to deter¬ 

mine the indices of domains in which the points (x, t) 
and (x,0) lie. Since 4>(x,f = 0) = 0, however, it suffices 
to know the number of domains s we shift to the left or 
to the right while we travel along the straight vertical 
segment S connecting the points (x,0) and (x,f). The 
number s is given by the difference of the numbers of 
rays intersecting the segment S = [(x,0), (x,f)] from the 
right and from the left, s = n + — n_. 

Clearly, at most one ray from each pair can intersect 
the vertical segment, and a ray of velocity v > 0 in¬ 
tersecting the segment from the left must be the right 
member of a pair originating from the spatial interval 
[(x — vt, 0), (x, 0)]. Since pairs are created uniformly at 
t = 0, the probability that the right going ray of a given 
pair with velocity v intersects the segment from the left is 
vt/L , where L is the size of the system. The probability 
that a given pair leads to such an intersection is 

r°° v t 

P = J 0 dv /L ^ v ^' ( 12 ) 


By symmetry, the probability of left intersections is the 
same. Since pairs are created uniformly, left and right in¬ 
tersections are independent Poissonian processes, and the 
probability of a pair configuration with n + (?z_) crossings 
from the right (left) is simply given by 

y{n + ,n_) = ——Q n + +n -e- 2 Q ( 13 ) 

n+!n_! 

with Q = tp / 0 °° dx vf(v) = (n±), and p = N/L the total 
density of pairs. 

Since at time t = 0 the soliton-antisoliton pairs shrink 
to single points, the domain l-\ of 4>(x, 0) lies between two 
pairs with probability 1 and l\ is even. If s = — ?z_ 

is even, then so is the domain I 2 = (li + s). In this 
case domains l\ and I 2 are separated by s/2 soliton- 
antisoliton pairs and have necessarily the same *1) values: 

— 4% = 4>(x,t) — 4>(x,0) = 0. If s is odd, however, 
then the Z| h domain lies at t = 0 + in a domain between 
the members of a pair in which 4>(x,t) = <f>(x, 0) ± 2tt//3 
with equal probability. Averaging over the two possibili¬ 
ties gives ( e l27ra //3 j re ~ l 2 ' KOt l ^)/2 = cos(27t a // 3 ). The final 
result for the expectation value is then 


G a 


00 

^2 p( n +’ n ~) 

n .|_ ,n_ =0 


/1 + (-i) n +- n 
V 2 


1 - (-i) n +- n - 
2 


cos(2to//3) 


(14) 


where Ga = ^e ra4> ( x,t l) i s the vacuum expectation 
value computed exactly in Ref. EH Carrying out the sum¬ 
mation yields then 

ia$(x,t)^ jG a = cos 2 (7ra//3) + sin 2 (7ra//3)e _t//T , (15) 
with the characteristic time r expressed as 

/*oo 

t~ 1 = 4p / d vvf{v) . (16) 

J 0 

The expectation value thus approaches an a-dependent 
constant exponentially fast, with the relaxation time in¬ 
dependent of a , that is, independent of the operator mea¬ 
sured. As a matter of fact, this time scale is, up to a 0{ 1) 
constant, given by the mean distance between the par¬ 
ticles p~ x divided by the average velocity, that is, the 
average time between two collisions. 

The somewhat surprising non-zero asymptotic value 
can be understood as follows. As explained earlier, do¬ 
mains conserve their ‘colors’ even upon collisions. There¬ 
fore, at any time, half of the domains must have phase 
4> = 0, while one quarter of them possess phases 4> = 
±2n//3, respectively. For times t r, solitons collide 
randomly and each individual domain exhibits a ran¬ 
dom Brownian motion. Therefore, at any given spa¬ 
tial point we find with probability 1/2 a phase 4> = 0 


while phases 4> = ±2tt/(3 occur with probabilities 1/4 
each. This immediately yields the asymptotic expecta¬ 
tion value /Q a —t cos 2 (to//?) 

Remarkably, expanding the result (14) for small t we 
recover the results of Ref. 1551 


-/ =l^sin 2 (TO/^)-+0(t 2 ), (17) 

Ga V 

obtained through the linked cluster expansion. There, 
however, higher order contributions being missing, a pure 
exponential decay to zero was assumed and an operator 
dependent decay rate, r” 1 = sin 2 (na / /3)t~ 1 , was de¬ 
fined. 

In the special case a = /3/2 + k tt, k G Z, the asymp¬ 
totic value is zero and a pure exponential decay is ob¬ 
tained. For these operators we can compute the expecta¬ 
tion value even without assuming f ++ (v) = f_ (v) = 0, 

i.e. in the presence of soliton-soliton and antisoliton- 
antisoliton pairs. In this case kinks crossing the segment 
S simply flip the sign of yielding 

(e^’*)/ 2 } = g N 2 e- t/T , ( 18 ) 

where now f(v) = YJmm 1 fmm'(v). We thus recover the 
exponential decay found in Ref. |55| (c.f. Eq. (|8|). No¬ 
tice that in Ref. 1551 the representative state approach was 
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based on a steady state computed in the leading order in 
the kink density but the series describing the time evolu¬ 
tion was resummed. It is thus a non-trivial fact that the 
semiclassical approach completely reproduces the result. 

For a = f3 + kn we get = Qp, so the operator 

does not evolve in time in the semiclassical approx¬ 
imation. This is probably due to the fact that it is the 
interaction energy density and the kinetic energy given 
by the derivatives are suppressed in the small quench 
limit. We note that the form factor calculatiorP^leads to 
the same trivial result. 


IV. RELAXATION OF CORRELATION 
FUNCTIONS: DERIVATION 

A. Case a = jH/2 

Before turning to the general case, let us compute the 
dynamical two-point function (|3j) of the operator 
with no restriction on the type of pairs. As we have seen 
in the previous section, this is a particularly simple case 
since 

e i/3/2[4>(x,t)-$(x',t')] _ = ( —l) n , (19) 

where n stands for the total number of trajectories in¬ 
tersecting the segment S = [(.T,f), (x',t')] connecting the 
two operator insertion points. This feature is the reason 
why we can allow general initial states with all possible 
kinds of kink pairs. 

Pairs of which both rays intersect the segment do not 
contribute. Let us compute the probability q that exactly 
one ray of a given pair crosses the segment, as shown in 
Figs [3bJ [3c] and [3f] Without loss of generality, we shall 
assume that x' > x. It will be useful to define the veloc¬ 
ities 


x — x 
t! T t 


x — x 
t' — t 


( 20 ) 


It is simple to check that if the velocity of the pair satisfies 
v < v then both rays of a pair can cross S. Then the real 
space domain where pairs with just one crossing ray can 
originate consists of two intervals, [x — vt, x + vt] and 
[x — vt', x + vt'] of total length 2 v(t' +1). In the opposite 
case, v > v, at most one ray of a pair can cross S. One of 
the two intervals in this case is [x — vt, x' — vt] for v < v s 
(Fig. 3b) or [x — vt, x' — vt] for v > v s (Fig. [3f|, and the 
other one is [x + vt, x' + vt'] (Fig. 3c). The lengths of the 
intervals are | (x' — x) — v(t' — t ) \ and |(a/ — x) + v(t' — t)\, 
where the modulus ensures that the expressions cover 
both the t! > t and t' < t cases. Thus the probability that 
one of the two rays of a pair with velocity v intersects the 
segment is 

2 v(t + t’) _. ]Ax — vAt\ + |Ax + vAt\ 

q v = Q(v-v)^- - J -+Q{v-v ) 1 -^- 1 

L/ Li 

( 21 ) 


with Q(x) the Heaviside function, Ax = x' — x, At = 
t' — t. Then the probability that only one ray of a given 
pair will cross is q = f Q d vf(v)q v , and the weight of a 
configuration having n crossings has Poisson statistics, 
p(n ) = Q n e~®, with Q expressed as 


Q(Ax', t, t') = Nq = 2p(t + t') f d vf(v)v 

Jo 

nV s /»00 

+ 2pAx / d vf[v) + 2p\t! — t\ / d vf(v)v (22) 

J V J V s 

and f(y) = f a b(v)- The correlator is then obtained 
by averaging |19|) over all values of n, yielding 

Cp /2 (Ax-,t,t')=C$ 2 (x'-x-,t'-t) , (23) 

with CfJ 2 (x' — x; t' — t) denoting the vacuum correlator. 

As we noted before, this expression should also de¬ 
scribe the correlation of the order parameter in the con¬ 
tinuum Ising field theory. Indeed, it agrees with the scal¬ 
ing limit of the exact resultPS for the transverse field Ising 
chain once an appropriate f(v) function is used in it. Nat¬ 
urally, it can also be obtained by taking the scaling limit 
of the semiclassical result for the Ising chain presented in 
Ref. 1451 It is important to note that it is not a priori ob¬ 
vious that the scaling limit of the post-quench behavior 
of the spin chain coincides with the post-quench behavior 
of the field theory, because a sudden quench can excite 
high energy states. In the semiclassical method it is how¬ 
ever almost automatic, as the derivations are essentially 
identical, and only the dispersion relation and the dis¬ 
tribution function differ. So, at least for quenches in the 
semiclassical regime, it appears that the field theory cap¬ 
tures correctly the non-equilibrium behavior of the spin 
chain. This was also observed in Ref. ED where it was 
found that the asymptotic time evolution of the order 
parameter after a mass quench within the paramagnetic 
phase of the Ising field theory agrees with the scaling 
limit of the time dependent magnetization of the Ising 
spin chain after quenching the magnetic field. 


B. General a 


Let us turn to the calculation of the time dependent 
correlation function © for general a. For the sake of 
simplicity, let us assume that x' > x, t' > t, since the 
other cases follow by symmetry. Unlike the previous sub¬ 
section, we now restrict the analysis to the case when 
there are only soliton-antisoliton and antisoliton-soliton 
pairs present with the same velocity distribution. As we 
mentioned, initial states with fixed $ values fall into this 
class. 

For a generic vertex operator (a arbitrary), the calcu¬ 
lation of the correlation function is similar to that of the 
expectation value |2|, but technically considerably more 
involved. According to Eq. (101 we need to determine the 
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(a) Avoiding pair (LL). 



(b) Left cut (L). 




(c) Right cut (R). (d) Double cut (D) 




FIG. 3: Representative examples of the classes of pairs 
defined in the main text. 


distribution of the number and orientation of domains 
between point (x, t) and (a/, t'). Depending on their rela¬ 
tion to the segment S = [(x, t), (x r , t')] we divide soliton- 
antisoliton pairs into six disjoint classes (representative 
examples are shown in Fig. [3): The first class contains 
avoiding pairs with both rays lying on the right (RR) or 
on the left (LL) of segment S (an LL pair is shown in Fig. 


3a). These pairs do not affect the correlation. We call left 


crossing (L) the pairs whose right going ray has a trajec¬ 
tory that crosses from the left and from below, while the 
left going ray of the same pair avoids the segment (see 
Fig. 3b). These pairs necessarily have velocities v < v s . 
We define in an analogous way right crossing pairs (R) 
which can have arbitrary velocity (see Fig. 3c). The class 
of double crossing pairs (D) includes pairs with both rays 
crossing S (see Fig. [3d|. We call pairs such that the seg¬ 
ment lies in between the two rays inclusions (I). Finally, 
a crucial role is played by those pairs whose right going 


ray crosses the segment “from above” (A) (see Fig. [3f| . 
This is only possible if the velocity of the pair is greater 
than v s . 

The probabilities that a given pair belongs to either of 
these classes are calculated based on the uniform spatial 
distribution of the pairs, in the same way as in the pre¬ 
vious sections. For example, the probability that a given 
pair belongs to class l L’ for t' > t is 

dv2vtf(v) + -j- J d v[x'- x-v(t'-t)]f(v) 

( 24 ) 

with f(v) = /_(_( v) + /_|-(u) = 2/-|_ (v). The probabili¬ 

ties for the other classes are listed in the Appendix. 

In a given configuration, the number of domains be¬ 
tween points (x,t) and ( x',t ') can be expressed in terms 
of the number of pairs in each class. Counting the do¬ 
mains from the left, the point ( x,t ) is in domain 

1 = 2n LL + n L + nj + 2 ha ( 25 ) 

while (x',t r ) is located in domain l 1 = 2nLL + 2.riL J rri R + 
2ri£> + nj + riA- Thus the number of kinks between the 
two points is 


s = l' — l = til + n R + 2no — ha ■ 


(26) 


Let us now focus on the average over the charges of these 
s kinks. We first remark that in this universal semiclas- 
sical limit kinks are strongly correlated: the charges of 
kinks 2 k — 1 and 2k always add up to 0. Therefore, if 
s is odd we have to average over the charge of a single 
kink without a pair, yielding ( e l27ra /^ e _l27r “/ /3 )/2 = 

cos(27ra//3). If, on the other hand, s is even, then there 
are two possibilities: either both l and l' are even, so 
that all s kinks form s/2 pairs and — <!>;= 0, or both 
l and V are odd, in which case we have s/2 — 1 pairs with 
zero total charge, and two unpaired kinks. Averaging over 
the charges of these two kinks yields cos 2 (27ra//3). We 
have to separate, however, the special case, s = 0: then 
$ii — $i = 0 irrespectively of the parity of l. Notice that 
the s = 0 configurations can be non-trivial due to the 
existence of the special cuts from above (A) that shift 
back the domain of (see the sign in Eq. (26) and 

Fig- [3fJ. 

Since only the parity of s and l matters, averaging over 
the velocities and positions of the pairs translates into a 
multiple sum over the numbers of the various types of 
pairs weighted by the probability of such a configuration: 


C a [x' - x;t,t') = ^2 


nR,n L ,riA,nD,ni 


_ Q n RQ nL Q n AQriD Qni ^—Qr—Ql — Qa — Qd — Qi 

n R \n L \n A \n D \ni \^ R L A d i 


( s , . , f 1 — (—l) s (2na\ l + (-l) s 

( 5 s ,o ■ 1 + (1 — Sgfi) ■ | - cos f —g— J 4 - 


1 + (- 1 )' ■ ! + 1 - (y 1 ) 1 cos 2 j 


cos 


(27) 
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where 


C a {x' — x; t, t') = C a (x' — x ; t, t')/C^\x' — x\t' — f) (28) 

is the correlation function normalized by the vacuum correlation, and Q^{Ax\t,t') = IVg^Aa;; t, f') with the proba¬ 
bilities q^ listed in the Appendix. 

The terms not multiplied by the Kronecker <5 Sj o can be written, collecting the sign factors and using some basic 
trigonometric identities, as 

Ci = cos 4 (na/P) + sin 4 (ira/P) e ~ 2 QR- 2 Q L ~ 2 Q A + sin 2 (na/P) cos 2 (ira/p) e^ 2 ® 1 (e~ 2 ® L + e~ 2 Q R ~ 2 ® A ). (29) 

The terms proportional to S St o can be dealt with by using the integral representation for the Kronecker delta, 5 Sj o = 
f* n and yield a contribution 


C 2 = 2sin 2 (™//3)cos 2 (7ra/^) U^-VQr. - P -^+i)Q,- 2 Ql 

J-k 2tt V 


(30) 


The total expression for the correlator normalized by its vacuum value is then expressed as 


C a (x' - x\t : t') = Ci + C 2 

= cos 4 (na/P) + sin 4 (na/P) e ~ 2 QR- 2 QL~ 2 QA s j n 2 (^ a //3 ) cos 2 (na/p) e~ 2Ql (e~ 2 ® L + e -2 ^^ -2 ^- 4 ) 


+ 2 sin 2 (tt a/P) cos 2 ( na/P ) e 


-( Ql+Qr+Qd+Qa ) ^tpQRe^+QDe^+QAe- 1 * ( p QLe i{t> _ „-Qie^-2Qj 


27T 




(31) 


We remark that expanding to second order in a and differentiating with respect to x and x' yields the topological 
charge density correlation function. 


V. RELAXATION OF CORRELATION 
FUNCTIONS: DISCUSSION OF THE RESULT 


We shall now analyze the general result in Eq. (31) and 
examine its behavior in various limits. First we note that 
setting a = P/2 + kn only the second term of Eq. (31) 
survives, and Eq. (31) simplifies to 


C/3/2 (^ 


x-t,t') = e - 2 Qfl- 2 Qi.- 2 Q^ ^ 


(32) 


and we recover the result (23) with 
defined earlier in Eq. (22). For a = 


Q — Qr + Ql + Qai 

P + kir, on the other 


hand, only the first term remains, and the semiclassical 
calculation yields Cp(x,x';t,t') = 1. 

Let us now analyze the local correlation function, the 
equal time correlation function, and the asymptotic dy¬ 
namic two-point functions for late times for general val¬ 
ues of a. 


A. Local correlation function 


The local correlation function is obtained by setting 
x —> x'. Using the special limits of the Q listed in Eq. 


(A8) of the Appendix, we obtain 


(70,(0; t, t') = cos 4 (na/P) + sin 4 (na/p) e At ^ T 

+ sin 2 ( 7 ra/P) cos 2 (ira/P) 


e“ t / T ( 1 + e~ At/T ) + 2e -A C( 2T )(i - e~ t/T )I 0 


, (33) 


where At = t' — t and r is the characteristic relaxation 


time defined in Eq. (16). 


r 


The local correlation function is computed for a veloc¬ 
ity distribution 


f{v) 



0 -v 2 /v 


2 

0 


(34) 






















FIG. 4: Local correlation function for different values of 
t as a function of t'/ t, where r is defined in Eq. (16). 
The decay is exponential for t = 0 and diffusive 

~ l /\/|t ~ A f° r all t > 0. 


and displayed in Fig. [4] This form of f(v) is motivated 
by the observation that the amplitude K(9) is an odd 
function of 9 around 6 = 

It is easy to check that C a (0\t,t) = 1 as At -> 0, a 
condition giving rise to the spike structures in Fig. [4] The 
most interesting feature of the result is its behavior for 
large time separation. For a finite and fixed t and large 
values of t! we find 


C a (0;t,t' —> oo) = cos 4 (nct/p) 


sin 2 (7ra/ /3) cos 2 (7ra/ P) 




2(1 


- P~ t / T \ 


yj TrAt/l 


(35) 


For any t > 0, the late time behavior is thus diffusive , 
similar to the behavior found within the semiclassical ap¬ 
proximation in thermal equilibrium 42 44 . As At —> oo, 
the local correlation function approaches a ^-dependent 
non-zero constant, equal to the product of expectation 
values given in Eq. (14) at times t and t' —¥ oo, as ex¬ 


pected. (The normalization factor given by the vacuum 
correlator also factorizes in this limit.) 

The origin of the diffusive contribution can be under¬ 
stood in terms of the simple picture of randomly diffusing 
domains, discussed in Sec. |III| as follows: for sufficiently 
long times l>r the point (0, t) is in a randomly selected 
domain l. This domain follows a diffusive motion, and re¬ 
mains at the point x = 0 with probability ~ l/V DAt, 
with D ~ l/{rp 2 ) the diffusion constant. With this prob¬ 
ability, the phase of the vertex functions at (0, t) and 
(0, £') is exactly the same, and the correlator gives 1. 


Interestingly, the diffusive term vanishes only for t = 0 
in which case we recover the exponential behavior found 
for the expectation value: 

C a { 0; 0 ,t r ) = cos 2 (ira/p) + sin 2 (na/P) e -t / T . (36) 

It is also interesting to note that for a = (3/2 + kir, 


C/ 3/ 2 (0 ;£,£') = e 


— p~ At / x 


(37) 


is independent of the time after the quench. This is true 
also for more general initial states having pairs with 
non-zero total charge, as can be seen from Eqs. ( 22|23 ) 
where the first two terms in Eq. (22) vanishes upon 
v s ,v —>• 0. This instant relaxation might be an artefact 
of the semiclassical approximation. However, the calcu¬ 
lation for a = P/2 is equivalent to that in the Ising 
field theory which can be mapped to free fermions, and 
in such systems similar behavior has already been ob¬ 
served. In Ref. [29( for a non-relativistic gas of bosons 
the exact density-density correlation function was ob¬ 
tained analytically after quenching the contact interac¬ 
tion strength from zero to infinity, exploiting the map¬ 
ping between the infinitely repulsive Bose gas (Tonks- 
Girardeau gas) and free fermions. The local correlation 
function was shown exactly to be time-independent. Sim¬ 
ilar behavior was found for relativistic free field theories 
in a different situation, after connecting two semi-infinite 
systems thermalized at different temperaturetPP Thus in 
the case of the Ising model, the time independence of the 
spin-spin local correlation function found in the semiclas¬ 
sical approach can turn out to be an exact result. 


B. Equal time correlation function 


Let us turn now to the relaxation of the equal time cor¬ 
relation function. The Q functions for t = t' are given in 
Eqs. (A10). Since Qa = 0, the (f -integral in Eq. (31) can 
easily be evaluated and the correlation function becomes 


C a (Ax-,t,t) = cos 4 (7ra/^)+sin 4 (na/P) e 4 p Ax + 4 Qo _|_ 
2 sin 2 (naff}) cos 2 (na/P) e Q n ~ 2 P Ax + e - */ T (l — e^^ D }) 


where 


<5d(Ax,£) = p f dv(Ax — 2vt)f{y). (39) 
Jo 

It can easily be shown that for Ax = 0 (Qd = 0) we 
recover again C a (0;t,t) = 1. 

The correlation function is plotted in Fig. [5] using the 
distribution in Eq. ( |34| . At any finite t > 0, as Ax —> oo 
it approaches exponentially a constant, 


C a (Ax —► oo; t, t) 


cos 2 (na/p) + sin 2 (na/P) e t ^ T 

(40) 
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after the quench. The connected part of the correlation 
function is thus exponential ~ e~ pAx with correlation 
length given by the density of kink pairs. 

A different result is obtained by taking the t oo 
limit first, where we obtain the asymptotic steady state 
correlation function, 

C a (Ax; t —► oo) = [cos 2 {na/p) + sin 2 ( 7 ra/p) e~ 2pAx ] 2 . 

( 11 ) 

The correlation length in the asymptotic state is thus 
£as = 1/(2 p) rather than 1/p, showing that the large 
separation and large time limits do not commute. 


FIG. 5: Equal time correlation funtion at different times 
as a function of x/d where d = l/(2p) is the mean 
interparticle spacing. 


This is the square of the expectation value (14) at time t, 
so the cluster decomposition property holds for all times 


C. Asymptotic steady state 


Finally, let us discuss the correlation function in the 
asymptotic steady state. Asymptotically, time transla¬ 
tional invariance is restored, and the probabilities given 
in Eqs. (All) become functions of the spatial and tempo¬ 
ral separations only. Consequently, the steady state cor¬ 
relation function can be expressed as 


{Ax; At) = lim C a {Ax;t,t') 

OO 

= cos 4 {na/p) + sin 4 {na/P) e~ 4pAx ~ A( ^ A + 2sin 2 {na/P) cos 2 {na/P) e^ 2pAx ~ 2QA I 0 (2 \JQa{Qa + 2pAx)\ . (42) 


with 

/»00 

( 3(4 (Ax, At) =P dv(vAt, — Ax)f(v) . (43) 

J V s 

Assuming a power law behavior f{v) = f$v k + 0(v k+l ) 
for small velocities, the asymptotic value is approached 
as ~ 1/P k+1 \ 

For infinite separation, Ax — > 00 , the cluster decompo¬ 
sition holds and we recover the square of the asymptotic 
value of the 1-point function, 

(Tf (Ax -a 00 ; At) -a g 2 a cos 4 ( 7 ra/P) = . 

(44) 

For large temporal separation, At /{At) pAx, Qa —> 

At /{At) — pAx, and to leading order 

(7“ s (Ax; At —> 00 ) = cos 4 (na/p) 

+ 2sin 2 {na/P) cos 2 {ira/P) ^ =, (45) 

\] 7tA£/t ’ 

that is, the diffusive behavior of the local correlation 
function is recovered. 


VI. SUMMARY AND OUTLOOK 

In this paper we developed a semiclassical approach 
to study the time evolution of correlation functions after 
small quenches in the repulsive regime of the sine-Gordon 
field theory. Assuming that the velocity of quasiparticles 
(kinks) is very small compared to the gap and using the 
asymptotic S-matrix, we were able to derive analytical 
approximations for the 1-point and general dynamical 2- 
point functions. Remarkably, this simple method — not 
restricted to integrable systems — allowed us to recover 
results for the 1-point functions obtained earlier by means 
of exact form factor expansions. Though these calcula¬ 
tions were performed in the repulsive regime, we expect 
our results to remain valid even in the attractive regime 
P < 1/ y/2 where breathers (bound states of solitons and 
antisolitons) are also present. Passing a breather does not 
change the value of < F, and in the low energy limit the 
scattering between the breathers and the kinks is purely 
transmissive. Thus breathers decouple from kinks in the 
semiclassical limit and should not influence the correla¬ 
tion functions dramatically. 

Our semiclassical calculations yield certain results that 
are somewhat surprising at first sight. We find, e.g., that 
the expectation value (e la4> ( x,t )) approaches exponen- 
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tially a constant value for generic cds (see Eq. ©)■ This 
is a consequence of the fact that the low-energy limit 
of the S-matrix is perfectly reflective, implying that the 
spatial sequence of the domains is conserved in time. Al¬ 
lowing also transmission in course of soliton-antisoliton 
scattering would break this pattern and lead to random 
domain sequences and a decay of the expectation value. 

We also obtained an analytic expression for 
the time evolution of the dynamic correlations 
^ e ia&(x,t) e -ia$(x' ,t')^ 0 f generic vertex operators 
and analyzed it in the various limits. We found that 
the cluster decomposition property holds during the 
non-equilibrium time evolution after the quench, at 
least for not very short times. The asymptotic value 
of the correlation functions is approached as a power 
law. The two-time correlations and the local correlation 
functions, in particular, show diffusive behavior. This 
can be understood heuristically based upon the picture 
of ‘magnetic’ domains performing a random walk. 

Both the diffusive correlations and the saturating ex¬ 
pectation values are consequences of the assumption of a 
perfectly reflective S-matrix. At any finite energy, how¬ 
ever, the transmission has a small but finite probability 
P t r ~ v 2 /c 2 , yielding a finite time scale, above which 
expectation values should decay to zero and the diffu¬ 


sive behavior is also expected to turn into an exponen¬ 
tial decay. In this sense our asymptotic results describe 
a sort of prerelaxation plateau which eventually decays 
at late times. We can estimate relatively simply the cor¬ 
responding time scale by considerations similar to those 
in Ref. 021 Within a time period T a given kink partic¬ 
ipates in T/t co \\ collisions with a collision time t co i oc r. 
So after time T ~ rc 2 /n 2 domains should change their 
color with a probability close to one and the reflective ap¬ 
proximation of the S-matrix should break down. These 
considerations point into a possible direction of improv¬ 
ing the present method by incorporating the leading non- 
reflective part of the S-matrix, or to use the full S-matrix. 
These explorations need, however, extensive numerical 
simulations and are beyond the scope of the present work. 
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APPENDIX 

In this Appendix we list the integral expressions for the probabilities that a given pair belongs to one of the classes 
defined in Sec. |IVB| These are calculated similarly to the probabilities computed in Sec. m and |IV A[ based on the 
uniform spatial distribution of pairs giving rise to rays that intersect or avoid the segment connecting the two operator 
insertion points in space and time. We work with the convention x' > x, t' > t and we use the notation 


v = 


x — X 
t f -\- t 


V s = 


X — X 

t' — t 


• Probability that a pair leads to a double intersection 

qD = j j J dv i x ' v ( t ' + *)]/(«)■ 

• Probability that a pair leads to an inclusion 

1 f Vs 1 f°° 

qi= -jr / d v[v(t' +t)- (x - x)]f(v) + — / Av2vtf(v). 


(Al) 


(A2) 


(A3) 


• Probability that a pair leads to one right intersection 

i r IT 3 . 1 f°° 

qR = — / dv2vt / f(v) + — / Av[x’— x + v(t' — + @(t' — t)— / d v[x'— x + v(t'— t)\f(y) 


(A4) 
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Probability that a pair leads to one left intersection 


i r i r- i r°° 

qL = Y dv2vtf{v) + — / dv[x'— x — v{t'— t)]f{v) + ®{t — t) — / dv[x — x — v{t — t)\f{v). 
L Jo 1 Ji L J Vs 

• Probability that a pair leads to one left intersection from above 

1 

q A = 0(t' - t)— J dv[x - x' + v{t' - t)\f(v). 

In the correlation functions these probabilities appear multiplied by the total number of pairs N, e.g. 

f-V 

Qd = Nq D =p dv[x' - x - v(t' + t)]f(v). 

Jo 

Below we list some limiting forms of these expressions. 

1. For the local correlation function x = x', so both v = v s = 0, and 


Qd = Ql = 0, Qr = Qa = 


if - t 
4 r 


Ql= 2 7' 


where 


/•OO 

r” 1 2 3 =4 p d vvf{v). 

Jo 

2. For the equal time correlation function t = t', the velocities become v = Ax/(21), v s = oo, and 

pv pOO j. 

Qd =P dv(Ax - 2vt)f(v), Qi = p dv(2vt - A x)f(v) = - pAx + Qd , 

Jo Jv 2 T 

nV f‘OC 

Ql = Qr = 2tp dvvf(v) + Axp dvf{v)=pAx-Q D , 
Jo Jv 


Qa = 0 ; 


(AS) 


(A6) 


(A7) 


(A8) 

(A9) 

(AlOa) 

(AlOb) 


where we used J 0 °° d vf(v) = 1. 

3. The asymptotic steady state is obtained by taking t, t' —> oo with At fixed, when we have v —> 0, and (for t' > t ) 


/•OO 

Qa = P dv(vAt - Ax)f(v), Qn-tO, Qj -> oo, 

J V s 

Qr = pAx + At/{A t), Ql = pAx - At/{A t) + Qa ■ 


(Alla) 

(Allb) 


The leading correction to these asymptotic expressions for large t comes from integrals of the type fj d vf(v){a + 
bvt ). Assuming a power law behavior f{v) = fov k + 0{v k+1 ) for small velocities, this gives a ~ l/t^ 1 ' 1 power 
law correction. 
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